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Abstract. Inflation creates perturbations for the large scale structures in the universe, but 
it also dilutes everything. Therefore it is pertinent that the end of inflation must explain how 
to excite the Standard Model dof along with the dark matter. In this paper we will briefly 
discuss the role of visible sector inflaton candidates which are embedded within the Minimal 
Supersymmetric Standard Model (MSSM) and discuss their merit on how well they match 
the current data from the Planck. Since the inflaton carries the Standard Model charges their 
decay naturally produces all the relevant dof with no dark/hidden sector radiation and no 
isocurvature fluctuations. We will first discuss a single supersymmetric flat direction model of 
inflation and demonstrate what parameter space is allowed by the Planck and the LHC. We 
will also consider where the perturbations are created by another light field which decays after 
inflation, known as a curvaton. The late decay of the curvaton can create observable non- 
Gaussianity. In the end we will discuss the role of a spectator field whose origin may not lie 
within the visible sector physics, but its sheer presence during inflation can still create all the 
perturbations responsible for the large scale structures including possible non-Gaussianity, 
while the inflaton is embedded within the visible sector which creates all the relevant matter 
including dark matter, but no dark radiation. 
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1 Introduction and motivation for a visible sector inflation 

The primordial inflation [1] is the simplest dynamical mechanism which explains the seed 
perturbations for the cosmic microwave background (CMB) radiation with almost Gaussian 
perturbations as suggested by the recent Planck data [2-4]. Since inflation dilutes everything 
other than stretching the initial vacuum fluctuations, after the end of inflation the coher- 
ent oscillations of the inflation must excite the Standard Model (SM) quarks and leptons 
at temperatures sufficiently high to realize SM baryons and dark matter in the current uni- 
verse [5, 6]. In principle, inflation could have occurred in many many phases [7], and perhaps 
even future-eternal [8, 9], but it must end in our Hubble patch with the right thermal history 
and three light neutrino species [2]. 

In principle inflaton whose potential drives inflation could be an arbitrary hidden sector 
field , whose properties can be constructed solely to match the observational data from 
Planck [3, 4]. However note that the CMB observables merely probe the gravitational aspect 
of the problem, it is not sensitive to the inflation's couplings to the SM matter and neither its 
origin. Typically hidden sector inflatons are SM gauge singlets 2 , whose mass and couplings 
can never be probed directly. In order to explain the universe filled with the SM quarks 
and leptons such an inflaton should primarily couple only to the SM sector [13], which is an 
ad-hoc assumption. A gauge singlet could in principle couple to other sectors, i.e. hidden 
or visible, there is no symmetry which can completely forbid their couplings to the hidden 
sector. 

Especially, string theory provides many viable SM gauge singlet inflaton candidates, for 
a review see [14]. Inflation is typically driven either by close string moduli or open string 

1 The sector of particle physics which does not carry the SM charges are typically called the hidden sector. 
Any beyond the SM (BSM) physics harbours many hidden sectors which may or may not couple to the SM 
sector or its minimal extensions. The latter sectors are usually called the visible sector as they comprise of SM 
fields or share the SM charges. Inflationary models can be constructed solely within a visible sector physics, 
for a review on inflation model buildings, see [5] . 

There could be more than one inflaton fields and perhaps even of the order of 0(1O 2 - 10 3 ) as in the case 
of assisted inflation [10]. Although there are some embeddings of such models of inflation within large SU(N) 
gauge theories [11], and in string theory [12], but it is highly unlikely that nature would prefer such a route 
since none of these fields can be embedded within a visible sector physics. 
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moduli. In either case they are SM gauge singlets - therefore it is not at all clear why and 
how such an inflaton would decay solely into the SM dof. Typically string compactification 
yields many moduli and hidden sectors [15]. A high scale inflation, i.e. scale which is higher 
than the mass of the moduli, could in principle excite all the moduli and dump all its entropy 
in the hidden sectors fields [16]. The reason for this is kinematical, the inflaton can decay 
into hidden sectors due to typically large branching ratio, i.e. there are more hidden sectors 
and only one visible SM sector. Furthermore, one might as well worry whether the inflaton 
could excite dark radiation, provided some of the dof become extremely light, such as in the 
case of string axions [17], or dark matter, such as in case of Kaluza-Klein dark matter [18]. 

In this respect, it is vital that the last phase of primordial inflation, i.e. last 50 — 60 e- 
foldings of inflation must end in a vacuum of BSM physics which can solely excite the relevant 
SM dof required for the success of Big Bang Nucleosynthesis (BBN), see for a review [19]. 
In this regard, Minimal Supersymmetric Standard Model (MSSM) [20] provides a perfect 
setup where all the matter content is known and can be probed at the LHC [21, 22] 3 . 
SUSY also helps inflation model building, since inflation needs a potential which remains 
sufficiently flat along which the slow-roll inflation can take place in order to generate the 
observed temperature anisotropy in the CMB. SUSY at any scale guarantees the flatness of 
such flat directions at a perturbative and a non-perturbative level (for a review see [24]), 
besides providing a falsifiable framework at low scales. Furthermore, the lightest SUSY 
particle can be absolutely stable under R-parity, and thus provides an ideal cold dark matter 
candidate [25]. 

The minimalistic realization would be to embed inflation, dark matter within MSSM 
which are all determined by the known SM couplings which provides credibility not only to 
particle physics but also to cosmological predictions. Our aim of this paper will be to show 
this within three paradigms for the inflationary cosmology - in all the cases inflation happens 
in the visible sector of MSSM. The origin of perturbations could be sourced from the visible 
sector physics or it might as well arise from the hidden sector, we will discuss the role of 
hidden sector here which might be responsible for creating mild non-Gaussianity. We will 
also discuss their merits in conjunction with the release of the Planck data along with the 
constraints arising from the LHC. 

2 Three paradigms within visible sector inflation 

One can envisage three realistic scenarios. A simple single field model of inflation and a 
scenario with multi-fields. In the latter case we can capture all the essence by mimicking 
just two fields - one which is inflaton and the other could be either curvaton [26-28], or a 
spectator field [29, 30] as the simplest examples. 

• A single field model of inflation: 

It is well known that a single field model of inflation with a canonical kinetic term 
will yield almost Gaussian perturbations. Of course, one can depart from the simplest 
assumptions to generate non-Gaissianity, such as sudden change in the potential, mod- 
ifying the initial vacuum from Bunch-Davis, or introducing non-canonical kinetic term, 
for a review on non-Gaussianity see [31]. All these have interesting consequences for the 

3 In principle the scale of SUSY could be wide ranging, if it is at low scales, such as 0(1) TeV, it can 
explain the hierarchy problem along with the possibility of directly detecting the SUSY partners of the SM 
fields at the LHC, for a recent review on naturalness, see [23]. 
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primordial non-Gaussianity, but most of them are severely constrained by the current 
observations [3] 4 . 

In this paper we will revisit the parameter space of a visible sector single field model 
of inflation embedded within MSSM with canonical kinetic term and with the Bunch- 
Davis initial vacuum condition [33-35] . In all these models inflation happens below the 
Planck scale and generate small non-Gaussianities. Furthermore, these models produce 
the right thermal history of the universe without any dark radiation. 

• Curvaton scenario: 

A light subdominant field during inflation can also seed the perturbations for the CMB. 
In the simplest scenarios it is assumed that the inflaton fluctuations are sub-dominant. 
The light field known as a curvaton [26-28] can slow roll after the end of inflation, and 
decays later on once the inflaton has completely decayed. While decaying the curvaton 
converts its initial isocurvature fluctuations into curvature perturbations. This conver- 
sion leads to a pure adiabatic fluctuations if the curvaton dominates while decaying, 
on the other hand if the curvaton decay products are sub-dominant compared to the 
energy density of the inflaton decay products, then there is a residual isocurvature 
fluctuations. Furthermore since the conversion itself is non- adiabatic, there is a genera- 
tion of non-Gaussian perturbations of the local configuration. In order not to generate 
residual isocurvature fluctuations, the inflaton decay products must thermalize with 
that of the curvaton decay products. A priori this is a non-trivial condition. The only 
way it could be satisfied provided both inflaton and curvaton can be embedded within 
the visible sector, i.e. MSSM, then this problem could be addressed amicably since 
both the fields would decay into the MSSM dof [36]. 

• Spectator scenario: 

This is a completely new paradigm where a light subdominant field like curvaton is 
present during inflation, but it decays into radiation much before the end of inflation [29, 
30]. The sheer presence of such a light field can create perturbations for the CMB, but 
since the field decays during inflation, its decay products need not be that of the SM or 
MSSM dof. In principle if inflation is occurring within a visible sector the perturbations 
can be seeded by the hidden sector field, which is advantageous for many theories of 
BSM including string theory. We will illustrate this for the first time with an example 
of inflation occurring within MSSM, while the spectator field is made up of arbitrary 
gauge singlet arising from the hidden sector physics. 

There could be two possibilities for the observed tensor to scalar ratio being negligible 
(r < 0.11 with 95% CL) [4]. The scale of inflation could be genuinely below the GUT 
scale, which is the case we will be considering in all the examples below, or the second 
option could be that the gravity is purely classical and so is the vacuum [37], while matter 
component is treated quantum mechanically; for a review on cosmological perturbation, 
see [38]. A linearized Einstein gravity has no source term, therefore for a classical gravity 
without any source for exciting gravity waves in a homogeneous and isotropic universe, the 
resultant primordial gravitational waves will be absolutely zero [37]. Any positive detection 
of primordial gravitational waves will indeed shed an important light on whether gravity 
should be treated classically or not. 

4 It is possible to obtain small bispectrum but large trispectrum in the CMB in the cyclic universe sce- 
nario [32]. 
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3 Inflection point potential for MSSM flat directions 



The MSSM provides nearly 300 gauge-invariant F-and D-flat directions [39, 40], which are 
all charged under the SM gauge group. Out of these flat directions, there are particularly 2 
-D-flat directions: udd and LLe, which carry the SM charges and can be the ideal inflaton 
candidates [33-35], where u, d correspond to the right handed squarks, L corresponds to 
the left handed slepton, and e corresponds to the right handed selectron. Both the inflaton 
candidates provide inflection point in their respective potentials where inflation can be driven 
for sufficiently large e-foldings of inflation to explain the current universe and explain the 
seed perturbations for the temperature anisotropy in the CMB [33-35]. 

Since both udd and LLe flat directions are lifted by higher order superpotential terms 
of the following form, which would provide non-vanishing A-term in the potential even at 
large VEVs, but below the cut-off scale: 

A $ 6 

(3.1) 

6M p 3 ' v ' 

where A ~ 0(1) 5 , and M p = 2.4 x 10 18 GeV is the reduced Planck mass. The scalar 
component of <E» superfield, denoted by tfi, is given by 6 

u+d+d L+L+e 

* = ^g-, ^ = ^—' (3 - 2) 

for the udd and LLe flat directions respectively. After minimizing the potential along the 
angular direction 9 ($ = <pe ie ), we can situate the real part of <f) by rotating it to the 
corresponding angles # m in- The scalar potential is then found to be [33, 34] 



V^)= 1 ml<p 2 -A^ + A 2 ^°, (3.3) 
v ; 2 6 M£ MS 



where and A are the soft breaking mass and the A-term respectively (A is a positive 
quantity since its phase is absorbed by a redefinition of 8 during the process). The masses 
for LLe and udd are given by: 

m~ + m~ + m~ 
ml = -± £ , (3.4) 



m~ + m~ + m~ 
d 

3 



ml = — / * . (3.5) 



The exact value of A is irrelevant for the CMB analysis, as it does not modify the CMB predictions. 
However it is possible to extract A by integrating out the heavy dof. In the case if the origin of these operators 
arise from either SU(5) or 50(10), then the typical value is of order A ~ O(10~ 2 ) for 5O(10) and A ~ 0(1) 
for 5(7(5), as shown in Ref. [41]. 

6 The representations for the flat directions are given by: uf = dj = ^tj<^, d~l = ~^% ( t > Here 

1 < a, ft, 7 < 3 are color indices, and 1 < i, j, k < 3 denote the quark families. The flatness constraints require 

that a / (5 / 7 and j 7^ k. Lf — \ ^ J ' ^^ = -73^o) ' ^ fe= -73^' wnere 1 < a , & < 2 are the weak 
isospin indices and 1 < i, j, k < 3 denote the lepton families. The flatness constraints require that a 7^ b and 
i / j 7= k. Note that the cosmological perturbations do not care which combination arises, as gravity couples 
universally. 
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Note that the masses are now VEV dependent, i.e. m 2 (</>). The inflationary perturbations 
will be able to constrain the inflaton mass only at the scale of inflation, i.e. (fio, while LHC 
will be able to constrain the masses at the LHC scale. However both the physical quantities 
are related to each other via RGE as we will discuss below. For 

A 2 

l + 4a 2 , (3.6) 



40m| 

v 2 



where a <1, there exists a point of inflection (</>o) in V(4>), where 



\xVw) 



1/4 

+ 0(a 2 ), (3.7) 

V"( ( f>o) = 0, (3.8) 

at which 

V{<j> ) = ^m 2 ^ 2 + 0(a 2 ) , (3.9) 

V'(<t> ) = 4a 2 m^ + 0{a A ) , (3.10) 

9 

TYl 

V'"(M = 32-^ + 0(a 2 ). (3.11) 

From now on we only keep the leading order terms in all expressions. Note that inflation 
occurs within an interval 7 

^ ° (3.12) 



60M2' 



in the vicinity of the point of inflection, within which the slow roll parameters e = (M 2 /2)(V' /VY 
and r\ = M 2 (V" /V) are smaller than 1. The Hubble expansion rate during inflation is given 
by 



H inf ~ ^_ 7/ . (3.13) 

In order to obtain the flat potential, it is crucial that the A((po)-term ought to be close 
to m^((j)Q) in the above potential Eq. (3.3). This can be obtained within two particular 
scenarios - (1) Gravity Mediation: in gravity-mediated SUSY breaking, the ^4-term and the 
soft SUSY breaking mass are of the same order of magnitude as the gravitino mass, i.e. 
rritf) ~ A ~ ""13/2 [44], and (2) Spilt SUSY: in Split SUSY scenario the scale of SUSY is 
high and sfermions are very heavy, the j4-term is typically protected by R-symmetry, see 
Refs. [45, 46], as a result the A-tevm. could be very small compared to the soft masses. 
However, if the Yukawa hierarchy arises from the Froggatt-Nielsen mechanism, then the 
A-term can be made as large as that of the soft mass, i.e. ~ A, as in the case of Ref. [47]. 



7 For a low scale inflation, setting the initial condition is always challenging. However in the case of a MSSM 
or string theory landscape where there are many false vacua at high and high scales, then it is conceivable that 
earlier phases of inflation could have occurred in those false vacua. This large vacuum energy could lift the 
fMSSM lat direction condensate either via quantum fluctuations [42], or via classical initial condition which 
happens at the level of background without any problem, see [43] . 
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Figure 1. (<f>o, m^) plane in which inflation is in agreement with the cosmological observations of the 
temperature anisotropy of the CMB fluctuations. The blue region shows the inflaton energy scale and 
inflaton mass which are compatible with the central value of the amplitude of the seed perturbations, 
P ( = 2.196 x 1CT 9 , and the la allowed range of spectral tilt n s = 0.9603 ± 0.0073 [2]. 



The above potential Eq. (3.3) has been studied extensively in Refs. [34, 48, 49]. The 
amplitude of density perturbations Sh arid the scalar spectral index n s are given by: 

4 « = ^=^^i^coB E ^], (3.14) 



and 



respectively, where 



n s = 1 - 4VA2 cot[A/boBEVA2], (3.15) 



A 2 = 900a 2 A/^ BE (^)\ (3.16) 



In the above, A/cobe is the number of e-foldings between the time when the observationally 
relevant perturbations are generated till the end of inflation and follows: A/"cobe — 66.9 + 
(l/4)ln(F(</)o)/M^) ~ 50. The running of the spectral tilt is negligible [34, 48, 49] within the 
current bound of the Planck observations [4] . The perturbations are due to single canonical 
field, therefore one would not expect large non-Gaussianity from this model. The observed 
non-Gaussianity parameter denoted by /nl < 1 is bounded by the slow roll parameters, see 
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Ref. [31], and is consistent with Planck [3]. The scale of inflation is low enough that one 
would not expect any observed tensor perturbations in any future CMB experiments 8 . 

Instant reheating and thermalization [52] occurs when a single MSSM flat direction is 
responsible for inflation. This is due to the gauge couplings of the inflaton to gauge/gaugino 
fields. Within 10 — 20 inflaton oscillations radiation-dominated universe prevails, as shown 
in Ref. [53] . The resultant reheat temperature at which all the MSSM dof are in thermal 
equilibrium (kinetic and chemical equilibrium) is given by [53] 

T rh ~ 2 x 10 s GeV. (3.17) 

Since the temperature of the universe is so high, it immediately thermalizes the LSP pro- 
vided it has gauge interactions. The LSP relic density is then given by the Standard (thermal) 
Freeze-out mechanism. In particular, if the neutralino is the LSP, its relic density is deter- 
mined by its annihilation and coannihilation rates [41, 54]. 

The advantage of realizing inflation in the visible sector is that it is possible to nail down 
the thermal history of the universe precisely. At temperatures below 10 — 100 GeV there will 
be no extra degrees of freedom in the thermal bath except that of the SM, therefore BBN 
can proceed without any trouble within low scale SUSY scenario. This reheat temperature 
is marginally compatible with the BBN bound for the gravitino mass 7713/2 > O(TeV). It is 
also sufficiently high that various mechanisms of baryogenesis may be invoked to generate 
the observed baryon asymmetry of the universe. 

In Fig. 1 we have explored a wide range of the inflaton mass, rriA,, where inflation 
can explain the observed temperature anisotropy in the CMB with the right amplitude, 
P c = 2.196 x 10" 9 , and the tilt in the power spectrum, n s = 0.9603 ± 0.0073 [2]. The 
observables and n s have been shown by blue region. We have restricted ourselves to VEV 
below the GUT scale. Within the current parameter range the model provides negligible 
running in the tilt which is well within the observed limit. We have allowed a wide range 
for m,/, and <j>Q (the inflection point) to show that inflation can indeed happen within SUSY 
from low scales to high scale SUSY breaking soft-masses. High scale soft masses could be 
made compatible within split-SUSY scenario [47]. 

Using renormalization group equations the mass of the inflaton can be evaluated at any 
energy scales, thus providing connection between physics at the very high energies in early 
universe and experimentally probed scales at LHC. For the udd flat direction RGE is [41, 54]: 

dm 2 , \f 2 9 



AA 1/16,,, 8,, 

where p, = po = 4>$ is the VEV at which inflation occurs. For LLe: 

dm 2 l /3 9 

P~ r/ = o -Mo go H M?g 2 

p dp 6tt 2 \2 22 10 iyi 



(3.18) 



„dA 1 /3, r 2 9, r 2 



(3.19) 



In this paper we will mostly consider this scenario. In order to obtain large observable tensor to scalar 
ratio, r, we will have to embed this model within N = 1 supergravity (SUGRA). This will modify the potential 
with a large vacuum energy density besides providing SUGRA corrections to mass and A-term [50, 51]. It 
has been shown that it is possible to obtain r ~ 0.05 for both inflaton flat directions: udd and LLe [51]. 
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Figure 2. (i^o,m^) plane showing the inflationary parameter space that may be ruled out in a likely 
case if SUSY is not found below 1 TeV. Green region denotes the exclusion if the inflaton is udd and 
yellow is for the LLe case. The blue band shows the {^>o,m^) values which are compatible with the 
central value of the amplitude of the seed perturbations, Pt — 2.196 x 10~ 9 , and the 3er allowed range 
of spectral tilt n s = 0.9603 ± 0.0219 [2]. 



where Mi, M2, M3 are U(l), SU(2) and SU(3) gaugino masses, which all, assuming SUSY 
models which obey universality conditions like constrained MSSM (CMSSM) [55] , equate to 
m 1 / 2 at the unification scale, and g\, 52 and 53 are the associated couplings. To solve these 
equations, one needs to take into account of the running of the gaugino masses and coupling 
constants which are given by, see [44]: 

(3( gi ) = ai gf P(^£)=0, (3-20) 

with a± = 11/167T 2 , «2 = 1/167T 2 and at = — 3/16-7T 2 . 

Within CMSSM one can try to constrain the inflaton mass for udd and LLe similar to 
the analysis of Ref . [54] . The current LHC searches for SUSY particles put a stringent limit 
on squarks and sleptons, see [56, 57], and as a result on the inflaton mass as shown in Fig. 
2. 9 



4 Inflation and curvaton both embedded within MSSM 

The curvaton scenario [26-28] is an alternative mechanism for creating perturbations. In this 
scenario, the density perturbations are sourced by the quantum fluctuations of a light scalar 

9 A very similar analysis could be carried out for the MSSM Higgs inflation: jj,H u Hd [58]. 
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field a, known as the curvaton, which makes a negligible contribution to the energy density 
during inflation and decays after the decay of the inflaton field (j). The advantage of the 
curvaton mechanism is that it can generate non-Gaussianity [26, 59] in the primordial density 
perturbations and also significant residual isocurvature perturbations, neither of which are 
possible in the usual single-field inflation models. Both signatures are now well constrained 
by the current Planck data [2, 3]. 

If the curvaton does not completely dominate the energy density at the time of its 
decay, the process of conversion of initial isocurvature perturbations into adiabatic curvature 
perturbations can enhance the local form of non-Gaussian fluctuations by 

/nl ~ - , for r < 1 , (4.1) 
r 

where 

r = — ^— (4.2) 

Pa + p-y 

is the curvaton's energy density ratio at the time the curvaton decays [26]. Here p 7 is the 
energy density of the radiation as the decay products of the inflaton. 

However, if either the curvaton or the inflaton belongs to a hidden sector of BSM, they 
may decay into other fields beyond the SM dof. There is no guarantee that the hidden and 
visible sector dof would reach thermal equilibrium before the BBN takes place. In this case, 
residual anti-correlated isocurvature perturbations are expected to be in conflict with the 
CMB data, which constrain them to be < 5% [4]. If the curvaton belongs to the visible 
sector but the inflaton does not, a value of r ~ 1 would avoid this conflict [60], but would 
render any non-Gaussianity undetectable. Note that if r ~ 1 the curvaton is solely responsible 
for exciting all the SM dof, so it must carry the SM charges in order to avoid dark radiation for 
instance [60, 61]. The curvaton scenario lends strong support to a visible sector dark matter 
such as neutralino in the case of the LSP, because either from the decay of the inflaton or 
from the curvaton, the neutralino would thermalize with the rest of the plasma soon after 
its decay, and its final abundance will be determined by its annihilation and co-annihilation 
rates. 

Keeping all these constraints in mind we need to embed both inflaton and curvaton 
within a visible sector of BSM physics where they both decay into the SM dof. Let us 
consider the case where the inflaton, (j), and the curvaton, a, both originate from different 
saddle point directions which are orthogonal to each other at least at the lowest orders in an 
effective field theory 10 . The total potential is 

V tot = V(<P) + U(a). (4.3) 

Let us first discuss the origin of the curvaton, which we take to be an i?-parity conserving 
-D-flat direction of the MSSM. For the purpose of illustration we consider that to be LLe, 
which is lifted by the non-renormalizable operator: 



10 In Ref. [36] we embedded inflation and curvaton both within visible sector, i.e. within MSSM, for the 
first time. We found that if the scale of inflation was higher than the effective mass of the curvaton then the 
spectral tilt would tend towards flat spectrum, i.e. n s ~ 1. In this section we will illustrate with the same 
potential how to obtain a spectral tilt close to n s = 0.9603 and the required local non-Gaussianity in the 
range of the Planck data [3, 4]. 



- 9 - 



where A is a non-renormalizable coupling induced by integrating out the heavy fields at the 
intermediate scale, M*, which could be close to the GUT scale, i.e. M* ~ Mqut- The scalar 
component of the £ superfield and its mass are given by: 

(L + L + e) m| + m| + m§ 
o = 1 , ml = -± ^ . (4.5) 



where at the lowest order the potential along the a direction is given by similar to Eq. (3.3) 



n. 



1 9, ,o AX <7 6 



0" 



10 



W = HM 2 -^ W + A^, (4.6) 

where A ~ m a ~ 0(1 — 10) TeV, are the soft SUSY-breaking terms. We will assume that 
the curvaton rolls on a saddle point of the potential, i.e. A = ^/40m a , so the saddle point 
lies at 

a =( ™° ) M*. (4.7) 

We now turn to the origin of V((j)) within the MSSM. Let us consider a flat-direction 
orthogonal to the curvaton. If the curvaton is LLe, the inflaton could be udd direction. In 
which case both inflaton and curvaton are embedded within MSSM. We take the inflaton 
direction to be squarks, typically they are expected to be heavier than the sleptons: 

u + d + d 

*-^7s— (48) 

Note that udd and LLe remain two independent directions for the entire range of VEVs. 
This flat direction will also be lifted by the non-renormalizable operators. However, at larger 
VEVs the potential energy density stored in the udd direction will be larger than that of 
LLe, so it would be lifted by higher order terms: 

W=T^^. (4.9) 



^ri 3m Mi 

m>2 



The potential at lowest order would be: 



6 5 6 8 (j) 11 

X2 m! + X3 m! + X4 m! + 



(4.10) 



where • • • contains the higher order terms. Note that the A m in Eq. (4.9) are all non- 
renormalizable couplings induced by integrating out the heavy fields at the intermediate 
scale. At energies below the cut-off scale these coefficients need not necessarily be of 0(1). 

Potentials like Eq. (4.10) were studied in Refs. [62, 63]. For A 2 < A 3 < A 4 < A n < 0(1), 
they provide a unique solution for which the first and second order derivatives of the potential 
vanish along both radial and angular direction in the complex plane: dV/dcf) = dV/d(jf = 
d 2 V/d(j) 2 = d 2 V/dcj)* 2 = (a saddle point condition). For the first three terms in Eq. (4.10), 
it is possible to show that this happens when 

55 

A^ = ^A 2 A 4 , (4.11) 

"Note that this curvaton potential has the same origin that of MSSM inflation discussed in Eq. (3.3). 
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Figure 3. A schematic timeline is shown for the curvaton scenario. The red and blue curves corre- 
spond to the energy densities of the infiaton <fr and curvaton a (or their decay products) respectively. 
In both cases the fiat directions decay into MSSM relativistic dof in less than one Hubble time. 

at the VEVs: <fr = (fro exp [wr/3, iir, £57r/3], where 



Concentrating on the real direction, the potential energy density stored in the infiaton sector 
is given by: 



where (fro <C M*. 

With the infiaton and curvaton potentials, the universe would then undergo the evo- 
lution as indicated in Fig. 3. The curvaton a only ends slow roll some time after inflation. 
In our case both <fr and a will decay within one Hubble time of the evolution. Both carry 
the SM gauge charges for which thermalization is similar to the case of instant preheating as 
shown in Refs. [52, 53]. As is also shown in Fig. 3, we will use subscripts "e", and "c" 
to indicate the Hubble exit of the relevant perturbations, the end of inflation, and the time 
when a decays, respectively. 

During inflation, both fields are slowly rolling. The curvaton a remains subdominant 
so inflation is totally determined by the infiaton (fr, where the e-folds from the Hubble exit of 
relevant modes to the end of inflation is iVj.. Since inflation is held near the saddle point, we 
can approximate the infiaton energy to be nearly constant, V{<fro)- The curvaton's motion 
can be solved by integrating out its slow roll equation of motion 




(4.12) 




(4.13) 




(4.14) 



where 



~ dN ~ 



da 



U'(a) 



(4.15) 



Hdt 



3H 2 " 
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After inflation, the inflaton decays into radiation, whose energy density then satisfies 

p^V(cf> )e- 4N \ (4.16) 

where N 2 is the number of e-folds of universe expansion after inflation, till the curvaton a 
decays. 

From the violation of the second order slow roll condition, i.e. r/ CTC = —1, we find the 
curvaton slowly rolls after inflation for the number of e-folds N 2 , which satisfies 

a c = a - Aae- 4N \ (4.17) 

where 

< 4I8 > 

characterizes the typical "width" of the slow roll region of the curvaton field near the saddle 
point. 

The slow roll equation of motion for the curvaton can also be integrated out after 
inflation, although now the universe is dominated by radiation. Similar to Eq. (4.14), here 
we have 

rcd i = r dN (4 . 19) 

<r Je 

Therefore Eq. (4.14) and Eq. (4.19) fully describe the motion of the curvaton, before it ends 
slow roll. Since, m a 3> H* after the end of slow roll, a decays into radiation and the universe 
then evolves adiabatically 12 . 

When the curvaton receives the quantum fluctuations 5a* at the Hubble exits, its initial 
perturbation will be maintained during inflation. Its perturbation evolves according to the 
perturbed Eq. (4.14), which is 

5a* 5a P 

^f-^ = °- (4 - 20) 

After inflation, the curvaton's perturbations are converted into the curvature perturba- 
tions when the curvaton decays 13 . This can be specified by the perturbation in the number 
of e-folds of the curvaton's slow roll after inflation, 5N2, and the perturbation in the curva- 
ton field at the end of slow roll, 5a c . Together they should comply the same end-of-slow-roll 
condition for the curvaton, i.e. 5r\ ac = 0, which gives 

— = -4<5jV 2 . (4.21) 
c c - c 

The curvaton perturbation, 5a, also evolves after the end of inflation, according to the 
perturbed Eq. (4.19), with the relation 

5a c 5a e 



a' 



5N 2 . (4.22) 



c 



Combining Eq. (4.20), Eq. (4.21) and Eq. (4.22), we are able to solve the perturbation in the 
number of e-folds before the curvaton ends slow roll 

5N 2 = (4.23) 



12 The saddle point keeps the effective mass near the saddle point vanishing, but the bare mass term of the 
curvaton is of the order of the soft SUSY breaking term, i.e. ~ 1 — 10 TeV. 



13 



In fact in our case it happens soon after the curvaton ends its slow roll condition, by virtue of m a ^> H t . 
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Figure 4. Parameter space is scanned for the curvaton model for M* = 5 x 10 16 GeV. The allowed 
band is painted in green, giving a positive local bispectrum within /nl = 2.7 ± 17.4, the latest Planck 
observational constraint for 3o\ The initial condition for the curvaton cr* and the coupling constant 
A have been picked (according to Eq. (4.29) and Eq. (4.31)) to always match the observed central 
values for the power spectrum Pq = 2.196 x 10 -9 and spectral index n s — 0.9603 [2]. 



After the curvaton ends slow roll, the a field instantly decays into relativistic dof within 
one Hubble time and we can ignore its evolution. The e-folds of the radiation dominated 
era, N$, can be written as a constant plus a quarter of the logarithmic of the total energy 
density, i.e. 

N 3 = J log(p 7 + U c ) + const. (4.24) 

When the initial perturbation 5a* is present, it also changes the energy density at the time 
the curvaton ends slow roll. The perturbation in it's energy density, in this case, would be 
small compared to that of the radiation. This is both because the curvaton is subdominant, 
and because its potential is relatively flat around the saddle point. Therefore, the major 
contribution to SN3 comes from the perturbation in the radiation energy density, which 
comes from perturbing Eq. (4.24) as (according to Eq. (4.16)) 

5N 3 = ^-5JV 2 = -(1 - r)5N 2 , (4.25) 
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where 

Uc „ U(a )_ e4N2 



P 7 + U c V{(j)Q) 

is the energy density ratio of the curvaton when it decays. 

After taking into account the energy density perturbation at the curvaton's end of slow 
roll, the universe then enters an adiabatic evolution, and no super Hubble perturbations will 
be generated. Therefore the total perturbation in the e-folds of expansion is 

SN = 8N 2 + 8N 3 = N a 5a* , (4.27) 

where 

"•-w.- (4 ' 28) 

Therefore the power spectrum of curvature perturbation becomes 

P( = N ^ = ( ■*""• )'. (429) 
The spectral index of the curvature perturbation can then be calculated as 

n. = - 2 ±°Z^A. (4.30) 
a In k Act 

From the above equations, we find that the observed spectral index n s constrains the initial 
condition of the curvaton by 

ft — 1 

(J*=a + -^-Act. (4.31) 

The local bispectrum can be calculated, according to Eq. (4.23) and Eq. (4.26), as 

^ = ^« J-^ = J^L = A~I. (4.32) 
J 6iV2 6N a da* 6N a 6a* 6r r y ' 

We can scan the parameter space in m a and V(4>o), and calculate the possible local 
bispectrum, as shown in Fig. 4. 

In Fig. 4, we require the coupling A < 1, and the curvaton remains subdominant (r < 
0.1). These two constraints narrow the allowed parameter space to the green band. The local 
bispectrum typically acquires /nl ~ O(20). As a specific example, we pick the inflaton as 
having the parameters A2 = 10 -9 and A4 = 10~ 3 , with an inflation energy scale relatively low, 
at V(<po) = (2.4 x 10 9 GeV) 4 , with negligible curvature perturbations. A low scale inflation is 
helpful in obtaining the right tilt in the power spectrum for the curvature perturbations. For 
the curvaton mass m a = 7.4 TeV, we can acquire the observed power spectrum of curvaton 
perturbation by taking the coupling constant A = 0.012. The energy density ratio when 
curvaton decays is r = 0.086. This gives the local bispectrum /nl = 8.3. 



5 Spectator mechanism with a visible sector inflation 

The spectator mechanism is a new mechanism, which has been proposed recently in Refs. [29, 
30]. The perturbations are created by a sub-dominant field which decays during inflation, 
known as a spectator. The spectator field cannot modify the dominant inflaton dynamics, 
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but it can leave its imprint in the cosmological perturbations. The decay of the spectator 
field can create non-Gaussianity due to the conversion of the entropy perturbations into the 
curvature perturbations. The spectator field decays completely into relativistic species during 
inflation, thus leaving no residual isocurvature perturbations. All the matter is created by 
the decay of the inflaton field after inflation, therefore it is important that the inflaton sector 
must be embedded within a well motivated visible sector. For all practical purposes the 
inflaton field's perturbations could be assumed to be sub dominant as compared to that of 
the spectator's. 

In order to create the right thermal history, we provide a simple example of embedding 
inflation within MSSMxC/(1)b-l gauge group, where the latter is also gauged. A simple 
-D-gauge invariant flat direction which can be the inflaton candidate in our case is given 
by [35, 64]: 

WDhNUuL, (5.1) 

where h is the Yukawa coupling and, N, H u , L are corresponding right handed neutrino, 
Higgs and slepton superfields. Note that the above superpotential can generate Dirac mass 
for the light neutrinos if the scale of U(1)b-l breaking is of order O(TeV) and the Yukawa 
is h ~ KT 11 - 1(T 12 . 

The inflaton field <j> corresponds to the superpotential: 



N + H u + L 
V3 



(5.2) 



So its potential can be written as 



where A is the trilinear A-term and the soft SUSY breaking mass term for the flat direction 
is given by: 

ml = f L - , (5.4) 

Note that for A = Auia, there exists a saddle point for which V'(0q) = V"(0q) = 0. The 
saddle point and the potential are given by: 

= VS^ = 6 x 10- m, (™**L) , (5.5) 

Here m u denotes the neutrino mass which is given by m v = h(H u ), with (H u ) ~ 174 GeV. 
For neutrino masses with a hierarchical pattern, the largest neutrino mass is m u ~ 0.05 eV in 
order to explain the atmospheric neutrino oscillations [65]. In our case we will be investigating 
a range of the inflaton masses which can accommodate the right handed sneutrino mass close 
to the low scale supersymmetry, and also yield the correct neutrino masses: 

~ 1 TeV (h < 10~ n ) to ~ 100 TeV (h < 10" 12 ) . (5.7) 

The above range of h also guarantees the curvature perturbation contribution by the inflaton 
is negligible. 
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I 1 1 1 t 

* Phase I c Phase II e 

Hubble exit cr ends slow roll End of inflation 

Figure 5. A schematic timeline for the spectator scenario is shown. Two phases during inflation 
have been shown after the relevant modes have left the Hubble patch. The energy densities of the 
inflaton, the curvaton, and the total are drawn in red, blue and green respectively. The dashed green 
curves show how a perturbation in the spectator field may affect the universe's evolution. 



After inflation the inflaton NH U L will start coherent oscillations and will dump all 
its energy into the light relativistic species of MSSM and the lightest of the right handed 
sneutrinos. In fact the lightest right handed sneutrino could be the dark matter candidate 
if it is the lightest SUSY particle. Dark matter analysis with the lightest sneutrino has been 
performed in Ref. [64]. Since U{1)b~l is gauged, it will lead to a quick thermalization of 
all the relativistic dof with a reheat temperature very similar to the analysis of Ref. [53] . 
The reheat temperature will be roughly given by T r h ~ 

(3O/^) 1 /4y(0 o )V4 „ 10 8 GeV for 

~ 1 TeV [53, 64]. 

Let us now imagine that the spectator field has a simple potential arising from some 
hidden sector physics. For the sake of illustration we consider this to have a flat potential 
with a hyperbolic tangent profile: 

U{a) = ^ fl + tanh — ) . (5.8) 



2 V <V 

where Uq and o~q are constant parameters whose values we will scan to show how this sim- 
ple potential can explain the amplitude of the perturbations and also the observable non- 
Gaussianity of the local form. Note that Uq <C V(4>o). 

A typical timeline for the spectator scenario in our case is summarized in Fig. 5. We 
will use the subscripts "*", "c" and "e" to indicate the respective slices as the relevant 
perturbations from the spectator field leave the Hubble patch during inflation, the spectator 
ends slow roll, and the end of inflation. 

We can solve the background evolution of the spectator field a, typically a ends slow 
roll when 

V,c = *\ ' = -1. (5.9) 
V{(f>o) 
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(c) The local bispectrum of curvature perturbations 
/nl gives a broad parameter space surviving from 
the Planck observational constraints [3]. 




(b) The running of spectral index lies inside the cur- 
rent observational bound dn 3 /dlnk = —0.0134 ± 
0.0270 (3a) [2]. 
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Figure 6. The cosmological parameters for the spectator model. The yellow shaded regions are 
excluded due to multiple constraints. The green bands lie with the parameter space that gives a 
spectral index inside n s = 0.9603 ± 0.219, and the local /nl = 2.7 ± 17.4 from the Planck observation 
(3<r) [3, 4]. The red contour lines are for the values of the respective parameters. Here we have taken 
the pivot scale c-folding =45. 



From this we can solve a c as 

„ / 4:M^U 
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for a c > a®. Here the subscript "c" indicates at the time when a ends slow roll. We can 
also solve the background motion of the spectator field a before it ends slow roll, since the 
inflaton is dominating with a constant energy density. By integrating the slow roll equation 
of motion for a, we obtain a, as a function of N, the remaining e-folds of inflation, as 

„ / 4M 2 [/ (iV-iV c + l) 

where N c is number of e-folds of inflation from a ends slow roll to the end of inflation. The 
second and third order slow roll parameters for a then simplify to 

MgU" 1 
Va= V(M = ~N-N c + V (5 ' 12) 
M^U'U'" 1 

6 = y(0 o ) 2 = {N-N c + iy ■ (5 ' 13) 

For the pivot scale N = N*, the spectral index n s , the local bispectrum /nl, and the local 
trispectrum (/nl are determined by r\ ai( in this case, giving the leading order terms, see for 
details [29, 30] 

n '- 1 = 2n " = - N.- 2 N e + l > (5 ' 14) 

/NL = _5 6T = 6(N.-N c + l)r' (5 ' 15) 

_ 25(2^, 25 

54r2 " 54(iV,-iV c + l)2r2- ^ ibj 



where 



{V{4> c ) + C/(a c )) F(0o) ' 

is the energy density ratio of the spectator a at its end-of-slow-roll boundary. 
The power spectrum of curvature perturbations then, becomes 

Hi ( U* V = (N* -N c + l) 2 rV 



(5.17) 



Pc 4n*\M*Ui) 3nmy Q ' (5 ' 18) 



Therefore to achieve the observed amplitude for P^, this requires do to take the value 

a 2 _ (N* -N c + l) 2 r 2 



(5.19) 



When the inflaton model is given, the parameters iV* and V$ are fixed. The spectral 
index n s , the local bispectrum /nl, the local trispectrum <7nl> and the relative value (Tq/H* 
then only depend on the r, the energy density ratio, and N c , the number of e-folds from 
the spectator ends slow roll to the end of inflation. We can then parametrically plot their 
dependences on iV c and r in Fig. 6. 

One can see that the model predicts the spectral tilt, (as shown in green shade which 
depicts the 2a range,) the negligible running of the spectral tilt and the local bispectrum in 
the range observed by the current Planck data [2-4]. We have also shown the value of g^i 
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in Fig. 6. Since a decays into radiation there is no residual isocurvature fluctuations, which 
matches the data perfectly well. 

Since the origin of a field belongs to the hidden sector, it could arise from billions of 
hidden sectors of string theory, all the inflationary models which claim a successful inflation- 
ary cosmology could potentially act like a spectator field. One of the common feature for the 
spectator field is the flat potential and this can be achieved in many string motivated models. 
Now the advantage is that these stringy origins need not have to explain the matter content 
of the universe. The latter could be obtained from the visible sector model of inflation. 

6 Nonlocal bispectra from curvaton and spectator mechanisms 

So far we have discussed the non-Gaussianity of local type or in the squeezed limit. In 
this section we briefly discuss how to obtain the equilateral and orthogonal types of the 
bispectrum. The bispectrum for any shape /NL(ki,k 2 ) is in general defined as [31] 

B(ki, k 2 , k 3 ) = ~/ NL (ki, k 2 )(P(k a )P(k 2 ) + P(k 2 )P(k 3 ) + P(k3)P(ki)). (6.1) 
5 

Here B and P indicate the strengths of the two-and three-point correlation functions of the 
gauge invariant curvature perturbation £. They are defined as 

(C(ki)C(k 2 )> = 8vr 3 P(k 1 )5 3 (k 1 + k 2 ), (6.2) 

(C(k!)C(k 2 )C(k 3 )) = 8vr 3 P(k 1 , k 2 , k 3 )5 3 (ki + k 2 + k 3 ), (6.3) 

In Eq. (6.3), the delta function indicates ki, k 2 , and k 3 form a closed momentum 
triangle, because of this multi-point correlation function expansion. Depending on the shape 
of the triangle, the bispectrum /Ni,(ki, k 2 ) can take different values, corresponding to different 
types of the bispectrum. The most considered and also best observationally constrained type 
of bispectrum, the "local" type which we have discussed above, corresponds to the squeezed 
limit k\ ~ fc 2 3> fc 3 . For any single field slow roll inflation, the local bispectrum of curvature 
perturbations is constrained to be small by the almost scale invariant spectrum [66]. 

In principle, every other shape of the momentum triangle corresponds to a unique type of 
bispectrum, and is considered "nonlocal". Among them, however, what raises people's most 
interests are the "equilateral" and "orthogonal" types. The equilateral type comes from 
taking the equilateral momentum triangle where k\ = k% = fc 3 , and the excess equilateral 
bispectrum (compared to the local one) typically can be generated by models with non- 
canonical kinetic terms [67, 68]. The orthogonal type is constructed in [68], to account 
for the bispectrum contribution that is orthogonal/uncorrelated to the local and equilateral 
shapes, for the general single field inflation. For this reason, the orthogonal bispectrum is in 
general a linear combination of various shapes, and does not originate from a single-shape 
definition. 

It has been shown in the past that the non-canonical inflatons can generate nonlocal 
bispectra with various patterns [67]. If we embed such a non-canonical field in the de Sitter 
universe, this non-canonical field can still contribute to the nonlocal bispectra as a curvaton 
or a spectator [69]. The inflaton field dominating the de Sitter universe can still come from 
the visible sector, as demonstrated in Section 4 and Section 5, to be responsible for the matter 
production. 

The generic nonlocal bispectrum for the curvature perturbation £ for any nonlocal 
shape, written as /ij£ oc ^(ki, k 2 ), can be estimated as follows. In both the curvaton and 
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the spectator cases, the universe evolves almost adiabatically till the boundary where the 
curvaton or spectator changes its equation of state significantly. Before this boundary, the 
gauge invariant perturbation of the curvaton or spectator, which is defined as 

C(T (k) = -m-H S -^-, (6.4) 

Pa 

does not evolve after the Hubble exit of the relevant modes. Here a is the curvaton or the 
spectator field we are concerned about, p a is its energy density, ip(k) is the scalar perturbation 
in the metric, and here the prime means derivative w.r.t the conformal time. We have omitted 
the time dependence. 

Here we can define the nonlocal bispectrum for the gauge invariant perturbation Co-(k), 
as /NL(5(ki, k 2 ), similarly with Eq. (6.1) to Eq. (6.3) 14 . In particular, we have 

(Cr(kl)Cr(k 2 )> = 87T 3 P ff (k 1 )(J 3 (k 1 + k 2 ), (6.5) 

<C(ki)G(k 2 )C(k3)> = 87r 3 P a (k 1 ,k 2 ,k 3 )<5 3 (k 1 +k 2 + k 3 ), (6.6) 
fl^ki.ka.ka) = ^/^ ) ) (k 1 ,k 2 )(P (T (k 1 )P CT (k 2 ) + P (T (k 2 )P CT (k 3 ) + P CT (k 3 )P (7 (k 1 )). (6.7) 

At the spectator or the curvaton boundary, the quantum fluctuations in a then transfer 
to the curvature perturbations. In the simplest setup where we can assume that the inflaton 
contribution to the curvature perturbation is negligible, the total curvature perturbation 
becomes 

C(k)=rC CT (k), (6.8) 

where r is the energy density ratio of a compared to the total energy density when it reaches 
the boundary, i.e. at the time when the perturbations are converted. 

After the boundary, the spectator is quickly redshifted away and the curvaton decays 
into radiation which has the same equation of state with the inflaton's decay products. The 
universe then becomes adiabatic and C(k) does not evolve afterwards. As the observation 
confirms, we can assume the local bispectrum, determined by Eq. (4.32) or Eq. (5.15), is 
much smaller than the nonlocal bispectrum we are interested in. Therefore the nonlocal 
bispectrum of the curvature perturbations can then be calculated, according to Eq. (6.8), as 



,(nioc),, , v 5P(ki,k 2 ,k 3 ) 

1 2j " 6(P(ki)P(k 2 ) + ...) 

1 Sfl^.ka.ks) 



-ki-k2 



r 6(P T (k 1 )P CT (k 2 ) + 



-k. 



1 



r/S(ki,k 2 ). (6.9) 



r 



'NL(o-) 



The above equation also works for both the equilateral and the orthogonal types of bispectra. 

The multi-point correlation functions has been well studied for a noncanonical slow 
roll field in a de Sitter universe, showing possible large equilateral and/or orthogonal types 



14 It should be noted that here /^'^'(ki, k^) is the nonlocal bispectrum of a in the de Sitter universe 
dominated by an inflaton field, as opposed to a's nonlocal bispectrum when itself serves as the inflaton. 
Therefore, typically the £ CT or /^'^'(ki, k2) of the general single field inflation, may be applied as those in 
this discussion, but only after a scaling relation (as a function of r). 
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Table 1. Benchmark points considered in this study of MSSM inflation, MSSM curvaton and a 
spectator scenario with MSSM x U(1)b-l flat direction inflaton. "/" means the model predic- 
tion satisfies the latest constraint by the Planck data. None of these models produce observable 
isocurvature perturbations. 



Planck Constraints (lu) 


MSSM inflation 


MSSM Curvaton 


Spectator 


Tensor-to-scalar ratio 
r < 0.11 (95% CL) [4] 


Negligible, / 


Negligible, / 


Negligible, / 


10 y P ? = 2.196™ [2] 


/ 


/ 


/ 


n s = 0.9603 ± 0.073 [2] 


/ 


/ 


/ 


dn s /d\nk = -0.0134 ± 0.0090 [4] 


< -0.002, / 


/ 


/ 


f local = 2.7 ±5.8 [3] 


<1, / 


Constrained, / 


/ 


/X U = -42±75 [3] 


<1, / 


Constrained, / 


/ 


/X = -25 ± 39 [3] 


<1, / 


Constrained, / 


/ 


Relativistic dof [2] 


only SM 


only SM 


only SM 



of bispectrum [67]. The Eq. (6.9) then gives a simple enhancing relation for the nonlocal 
bispectrum of curvature perturbations, showing that the same mechanism also works in the 
curvaton or the spectator setup, in order to generate large nonlocal bispectra. 

Typically non-canonical kinetic terms such those as appearing in the DBI-innation occur 
in the hidden sector rather than the visible sector particle physics. It is simple to hide the 
uncertainties of the hidden sector physics in the spectator scenario as compared to the visible 
sector curvaton. One can imagine a spectator field with a non-canonical kinetic term slowly 
rolling its potential and decaying well before the end of visible sector inflation. The decay 
products of the spectator field would be diluted away anyway during the remaining inflation, 
rendering the universe solely with the visible sector inflaton which produces all the relevant 
matter for the BBN and without any trace of dark radiation. 

7 Summary 

In this paper we have emphasized only the visible sector models of inflation, curvaton and 
spectator scenarios in light of Planck data [2-4]. There are plethora of hidden sector gauge 
singlet models of inflation and curvaton with ad- hoc couplings and mass parameters, but 
these models are unable to explain why such a gauge singlet inflaton or curvaton would 
decay solely into the visible sector dof. We summarize our findings in Table- 1. Although we 
have restricted ourselves to the local form of /nl but one advantage of a spectator scenario 
is that it can generate nonlocal forms of /nl compatible with the Planck data. 

The spectator mechanism is amazing since it can fit all the observed parameters, 
n s , dn s /d\n k, f^ al , f§s^ h , /nl™' °^ ^ ne CMB without generating any isocurvature per- 
turbations, besides providing all the SM dof, which means no dark radiation by virtue of the 
inflaton being within the visible sector physics. 
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